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Abstract 

For one-matrix models with polynomial potentials, the explicit relationship between the partition function 
and the isomonodromic tau function for the 2x2 polynomial differential systems satisfied by the associated 
orthogonal polynomials is derived. 

1 Introduction 

It is well-known |34| that the partition functions for matrix models, both for finite n x n random 

matrices and suitable n — > 00 limits, provide r functions for the KP hierarchy. It is also known that such 
partition functions satisfy certain constraints amounting to the invariance of the corresponding elements 
in the associated infinite Grassmannian under a set of positive Virasoro generators. This generally implies 
that there exist differential operators in the spectral parameter variable which annihilate the associated 
Baker-Akhiezer functions. The compatibility of the differential equations in the spectral parameter and 
in the deformation parameters representing KP-flows implies that the monodromy of the operator in the 
spectral parameter is invariant under such deformations. In the large n limit this was considered in p3], 
while for finite n it was studied in |2l] . Associated to the analyticity properties of the Baker-Akhiezer 
functions is an isomonodromic tau function [ p2| , defined similarly to the KP tau function. 

In practice, the exact relationship between such "isomonodromic" tau functions and the KP tau 
function is not very clearly understood. There are special cases, however (see e.g. [Q, |[ ||[ [y| [ij], 
|2l| , |22] , |2(| ) where these two quantities can either be shown to be identical, or at least explicitly related. 
In the computation of spacing distributions for random matrices, a class of isomonodromic tau functions 
occurs, given by Fredholm determinants of "integrable" integral operators |)|, [l7], [l8[ |28, 32, 33| supported 



on a union of intervals, with the deformation parameters taken as the endpoints of the intervals rather 
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than as KP flow parameters. The precise relation between such deformations and the KP flows is not yet 
clearly understood. (See however, jl], [|, |l6[ [54| for some indications of the links between the two.) 

It is our purpose here to analyze in detail, for a generalization of the simplest version of the Hermitian 
1-matrix model, with measures that are exponentials of polynomial trace invariants, exactly what this re- 
lationship is. The main idea is to consider, for finite n (the size of the random matrix), the isomonodromic 
deformation systems satisfied by the associated sequence of orthogonal polynomials, and to compute an 



explicit formula (Thm. 2.4) relating the associated isomonodromic tau function to the partition function 
of the generalized matrix model. In the process, we recall (and generalize) some of the standard results 
concerning the relation of isomonodromic deformations to random matrices and the occurrence of Stokes 
phenomena in solutions to the associated systems of differential equations with polynomial coefficients. 

2 Orthogonal polynomials for semiclassical functionals 

We consider the orthogonal polynomials for a moment functional 

C^r ■ C[x] i — ► C , (2-la) 
Cxv{x l ) ■= [ dxe~ v{x) x l , (2-lb) 

where the integration over xT represents a suitable linear combination of integrals over contours to be 
explained in detail below. The function V(x) is called the potential and it will be taken here to be a 
polynomial of degree d+1, 

d+i K 

V(x) := £ u K — , (2-2) 

K=l 

in which the coefficients {ttjf} are viewed as deformation parameters. (The constant term is suppressed 
since it only contributes to the overall normalization.) 

This is a generalized setting which reduces to the ordinary orthogonal polynomials when the measure 
is positive and the contour of integration is on the real axis. We do not assume that V{x) is of even 



degree or that it is real; the moment functional (2-1) is in general complex- valued. In this general setting 



the orthogonal polynomials (if they exist) are defined as follows. 

Definition 2.1 The (monic) orthogonal polynomials {p n {x)} n eN are defined by the requirements 

£>cr(Pn(%)Pm(v)) = h n$nm , K E C X (2-3a) 

Pn (x) =x n + ... . (2-3b) 

It is easily seen |5lJ that such polynomials exist if and only if the Hankcl determinants formed from the 
moments do not vanish 

A n (xT, V) := det {C v {x^)) + , Vn e N . (2-4) 



(If they do exist, eqs. ( 2-3a , |2-3b| ) determine them uniquely.) 



We now return to the definition of the cycle xT . The convergence of the integral defining the moment 
functional C implies that the contour of integration must extend to infinity in such a way that the real 
part of V{x) diverges to +oo. For the purpose of the description of xT let us introduce the sectors 

^arg^+0 fc = Q 
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Figure 1: An example of the contours Tk in the case d = 7. 



These sectors are defined in such a way that $l(V(x)) — > +00 in the even numbered sectors. Therefore a 
contour should come from infinity in an even sector and return to infinity in a different even sector. Such 
functionals have been studied in a similar setting in ]23|, Q . 

We define the contours Tk, k € Zd as the ones coming from infinity in the sector J>2fc-2 and returning 
to infinity in the sector S2k (see Figure [j]). It follows from the Cauchy theorem that the corresponding 
d + 1 functionals Cr k satisfy 

d 
k=0 

but it can be shown that any d of them are linearly independent f23|| . In general we may consider any 
linear combination of such moment functionals corresponding to a "one-cycle" and the corresponding 
deformation theory under changes of complex generalized measure. The quantities xT and C^r arc 
understood to mean: 

d 



xT := ^ xfclfc , (2-7a) 
fc=i 

d d„ „ 

£ *r(0 :=5> fe £r fc (-) == $>* dxe' v ^(-)=: / dxe- y ^(-) , (2-7b) 



fc=i fe=i 



where {>ffe € C}fe = i ! ... ( j are any fixed set of constants (not all vanishing). (A similar setting, in which the 
contours Tk play the roles of boundaries between Stokes sectors, was considered, e.g., in the earlier works 

§0 §113. 

For a fixed potential V(x) the Hankel determinants A n (xT, V) are homogeneous polynomials in h of 
degree n. The nonvanishing of this denumerable set of polynomials determines a set with full (Lebesgue) 
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measure in the >c space and in this sense the orthogonal polynomials exist for generic >r. (An interesting 
problem, which will not be addressed here, is whether the set of x's for which the orthogonal polynomials 
exist contains an open set, thus allowing for arbitrary infinitesimal deformations.) 

If the contour of integration is not the real axis this setting does not derive in general from a matrix 
model. Nevertheless, we introduce what we call the 'partition function, in analogy with the Hermitian 
matrix model case pi, |3lJ , through the following formula 



Z„:= [dxf- Ux n \(xfe-^U V ^) = n \A n (x,V) , A(x) := J](a;, - Xj ) 



(2-8) 



It is well known that Z n is proportional to the product of all normalization factors for the monic orthogonal 
polynomials in Dcf. 



2.1 



n-1 

n! J] h, 



Instead of the p„'s, it is convenient to introduce the following quasi-polynomial functions: 

$ n {x) = ^=p n {x)<r* v W , 

Vi 



which are orthonormal, in the sense that 



dx lj)n{x)i> m {x) = Sum ■ 

r 

These (quasi-)polynomials satisfy a three-term recursion relation [ Ml : 

oo 

X1p n {x) = ^ QnrnMhn = Jn+l^n+1 + Pn^Pn + 

where 



m— 



7ti — Qn,n-1 



h n -l 



^ Vn S N, and /?„ 



(2-9) 



(2-10) 



(2-11) 



(2-12) 



(2-13) 



The semi-infinite matrix Q is symmetric (Q n m = Qmn), and has its nonzero entries along a diagonal 
band of width 3 {Q mn = if |m — n| > 1); i.e., it is a symmetric Jacobi matrix: 



Q = 



(Po 7i \ 

7i Pi 72 

72 73 

V ; o ■•. ■. -J 



(2-14) 



These recursion relations may equivalently be expressed as a 2-component vector recursion: 



where 



tp n+1 (x) 



RJx) 



Rn(x) 



1pn~l{x) 
4>n(x) 



1 

~f n X — n 

7w+l 7n + l 



(2-15) 



(2-16) 
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We can also express the derivatives of the ^ n 's in the basis {ipm} m = ^ 

oo 

l/>' n ( X ) = Pnm1pm{x) , (2-17) 

m=0 

where P has non-zero entries in a band of width 2d+ 1 centered on the principal diagonal. An integration 

(2-18) 



by parts using eq. (2-11) shows that the semi-infinite matrix P is antisymmetric. From 

p' n = n Pn . 1 +0(x n - 2 ) , 
one obtains: 

Pnm + \ (V'(Q)) nm = if m>n (2-19a) 



1 



P n ,n-1 + -X (V'(Q))n,n-l = ~ > 



and hence 



P 



In 



-~(V'(Q) + -V'(Q)-) 



(2-19b) 



(2-20) 



where for any semi-infinite matrix A, A + (resp. A_), denotes the upper (resp. lower) triangular part 
of A. Eq. ( 2-20D also implies the following equations, referred to sometimes as "string equations" 

[Q,P] = 1, = V'(Q) nn , — = V'(Q) n , n -i . (2-21) 

In 



2.1 Differential system 



It follows from (2-20) that the sum in (2-17) contains only a finite number of terms and can be written: 



1 d 

^) = - « £ ( V '(Q)n,n + ki>n+k - V'(Q) r 



l z4>n 



(2-22) 



fc=i 



Using eq. (2-12) recursively for any < k < d we can write tp n +k and V'n-i-fc as a unique linear 

combination of ij} n and tp n —i with coefficients that are polynomials in x of degree < k. Eq. 2-22 can thus 
be rewritten as a 2 x 2 differential system: 

d 



da; 



1pn(x) 



V n [x) 



1pn{%) 



n > 1 



(2-23) 



where T> n (x) is a 2 x 2 matrix with polynomial coefficients of degree at most d = deg V . In fact we can 
choose the sequence of 2?„'s to satisfy the recursive gauge transformation relations 

1 +d x R n R n - 1 . (2-24) 



— RnD n R 



n^n^^n 



(See also || for similar statements for the case of 2-matrix models.) This relation determines the sequence 
of T>nS uniquely ji]. |[ |l^, ^0|. |T). It follo ws in particular that these transformations preserve the 
(generalized) monodromy of the linear systems (2-23) for all n. 

Theorem 2.1 The matrix T> n (x) is given by the formula: 

' V'{Q)-V'(x) \ ( V'(Q)-V'( X ) 

1 fl 1 V Q~* ) „-i,„_i V 0-» 



V n (x) = -V'{x) 



-y'(x) 



1 o 

-1 

1 

-1 



+ 7« 



Q-x 

V'{Q)-V'(x) 
Q-x 

V'{Q)-V'(x) 



V'(Q)-V'(x) 
Q-x 



n—l,n 



V'(Q)-V'(x) 
Q-x 

V'(Q)-V'(x) 
Q-x 



~ln 





(- 



Q-x 



n— 1,7?.— 1 

71 ,71— 1 



(2-25a) 
(2-25b) 
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Proofs of this result may be found in j|, |L [l9| and earlier in plj for the case of even potentials. In 
the Appendix we give an alternative proof which is based upon the deformation formulas obtained in 
Theorem |2.2| below. 



-1 

1 

Darboux formula written here for the quasi-polynomials 



Remark 2.1 The matrix 7„ 



appearing in the first line of eq. ( 2-25b| ) arises naturally in the ChristofFel 
lasi-polynomials 

(x - x') ^2 4'j( X )' l Pj( X ') = 7n(lpn(x)lj} n -l{x') - ll> n -l(x)lj} n (x')\ = 7„ fa/>n-l (x) , ^ n (x)] 



3=0 



"0 -1" 




tp n -i(x') 


1 




Tpn(x') 



(2-26) 



Equation ( [2-24 ) is the starting point that allows us to explicitly compute the spectral curve of the 
matrix T> n (x) in terms of the deformation parameters {uk}k=i...<i and the partition function Z n . 

In addition to eq. (2-23) we may consider the effect of deformations of the {uk} parameters on the 
orthogonal polynomials. It is well-known that these give rise to isospectral deformations of the Jacobi 
matrix Q (Toda hierarchy equations Jl5|, |27j ) 



d 



du K 

d 
du K 

U K : 



Q = [U K ,Q] , 



1 

K 



{Q K + - Q K -) 



(2-27a) 
(2-27b) 

(2-27c) 



Combining the recursion relations (2-12) with equations (2-27a, 2-27b, 2-27c) one obtains deformation 
equations for the successive pairs of orthogonal quasi-polynomials 



d 



du 



K 



U K ,n{x) 



1pn{x) 

ip n+ i(x) 



(2-28) 



The matrices Uk,u (x) can be explicitly computed and chosen to have polynomial entries of degree < K 
as given in the following theorem, whose proof is given in the appendix. (Similar formulas for the case of 
even potentials were given in pjj . The main simplification in that case is that the coefficients j3 n vanish 
for all n's.) 



Theorem 2.2 The expression for Uk.u(x) may be taken as 



U K , n {x) 



1 

2K 



,.K 



in— l,n—l 



QK „k 

^Cnn 



K 



( * K -Q K \ 

\ *~Q )n,n-l 



x-Q 



( x K -Q K \ 

V x-Q Jn-l,n-l 

V ^ )n.n-l 



Remark 2.2 It follows as a corollary to Thm. 



2.2 



that: 



du K 



ln7n 



(2-29) 



(2-30) 



It is important to note that eqs. ( 2-15 ), ( 2-23 ) and ( 2-28 ) do not just hold for the sequence {ipn} of 
orthogonal quasi-polynomials but define an overdetermined system of differential-difference-deformation 
equations that are actually compatible in the Frobenius sense. Hence there exists, for all n, a basis of 



6 



simultaneous solutions of these equations. This system of solutions can be written explicitly by adding 
to the ipn's a solution of the second type as follows 



4> n -l(x) 1pn-l(x) 



i) n (x) = e? v(x) / da; ■ 

J xT 



(2-31a) 
(2-31b) 



It is an easy exercise to show that these ip n 's satisfy the same recursion relations as the ip n 's for n > d and 
hence satis fy th e same ODE (2-23). It is also straightforward to verify that they satisfy the deformation 
equations ( [2-28| ) . The fact that these equations are compatible implies that the generalized monodromy 
data (i.e. the Stokes matrices) for the operator d x —T> n (x) are independent of the deformation parameters 
{uk} and the integer n. 

We now turn to the main object of this paper: to relate the partition function (2-8) with the isomon- 
odromic tau function (whose definition is given below). The first step is to connect the spectral curve of 
the matrix T> n (x) with the partition function. 



Theorem 2.3 The spectral curve is given by 



= &ct{ y \-V n (x))=y 2 -\Tr(V n 2 {x))=y 2 -\{V\x)) 2 +^ 



3=0 



V\Q)-V'{x) 



Q - x 



d-K-l 



J=l 



K=0 



J=l 



\{V'{x)f + n £ u J+1 x J - 1 + g x K Uj + K + 2J ^rj ln ( Z ") 



(2-32a) 



(2-32b) 



Before proving this result, we remark that the RHS contains matrix elements of Q with all indices down 
to n = while, as it stands, formula ( 2-25b| ) for V n (x) contains only matrix elements of Q with indices 
"around" n, i.e., at distances less than d + 1 from Q nn . 

Proof of Thm. 2.3. First, notice that T> n (x) is traceless (by formula ( [2-25b )). Using eq. ( [2-24 ) and 
the cyclicity of traces of products one finds that 



Tr (V n+1 {x) 2 ^j = Tr (v n (x) 2 ) + 2Tr (V n (x)R n - 1 (x)R' n (x)) + Tr ^(R'^Rn^ix))' 
Now note that 



R' n {x)R n ~\x) = 





1 o 

In 



- + 





so the third trace term vanishes and the second one gives 



Therefore 



Tr \V n+l {xY) =Tx[V n {xY) -2 



n—i 

Tr(v n (x) 2 ) =Tr(v 1 (x) 2 ) - 2 £ 

3=1 



V'(Q) - V'{x) 



V'(Q) - V'(x) 



(2-33) 
(2-34) 

(2-35) 
(2-36) 



j j 



Note that T>q is not defined. We therefore need to compute the trace of the square of the first matrix T>i. 
To this end we introduce the notation 



W(x) := 



v'(Q) - YM 

Q - x 



(2-37) 
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Now use the string equation (2-21) to compute 

-V'(x) = [V'(Q) - V'(x)] 00 = [W(x)(Q - x)] 00 = (A) - x)W o(x) + jiW 10 (x) (2-38a) 



- = [V'(Q) - V'(x)} 10 = [W(x)(Q - x)] 10 = jiWnix) + (A, - x)W w (x) 



7i 



(2-38b) 



Multiplying eq. ( [2-38a ) by jiWip , eq. ( |2-38b ) by 71 Woo and subtracting, we obtain (since the diagonal 
elements of V'(Q) vanish by eq. ( 2-21 )) 



W 00 (x) = yi 2 W n (x)W 00 {x) - -fiW w (x) V'(x) + 7l iy M (x) 



Computing the trace of the square of T>\(x), and using eq. ( 2-39 ) we obtain 



Tr (V 1 (x) 2 ) = -V'(x) 2 -2W Q0 (x) . 



Therefore the full formula for the trace of the square is 



n — 1 

Tr(v n (xf) =-V'(xf-2j2 

j=0 



V'(Q) - V'(x) 
Q — x 



j j 



(2-39) 



(2-40) 



(2-41) 



Remark 2.3 Notice that, had we interpreted formula (2-251) for n = to mean that 70 = 0, we would have 
obtained exactly the same result. 



We now expand the polynomial 
V'(Q) - V'(x) 



d-2 



d-K-1 



j=i 



j=i 



a, J+1 x J 1 l+J2 xK £ u.J+K+2Q J 

K=0 

and use the fact (proved in the appendix, eq. ( |A-5| )) that 

Kd UK ln(Z„) , 



(2-42) 



3=0 



(2-43) 



to complete the proof of the theorem. Q.E.D. 



Theorem 2.3 relates the spectral curve of the matrix T> n {x) to the partition function Z n . The latter is 
known to be a r- function for the KP hierarchy. (A short proof of this fact is given in Appendix [b| together 
with an interpretation of Z n as a Fredholm determinant.) On the other hand, a linear system of the 
form (2-23) with compatible deformation equations (2-15), (2-2S), allows one to define an isomonodromic 
tau function p2[ . We recall here briefly, in a way adapted to the case at hand, how to define such a 
tau-function. 

The ODE 

^(x) = V n {x)*{x) (2-44) 
ax 

has just one irregular singularity at x = 00 and hence has Stokes sectors, which are just those defined in 
eq. (|2-5|). Since the leading term of the singularity in T> n (x) is semisimple with distinct eigenvalues 



V n {x) 



--±-x d a 3 + 0(x d ^) 



(2-45) 



we can write in each Stokes sector Sk a formal asymptotic solution of the form 

V n (x)~C n Y n (x)e T ^ , (2-46a) 

OO 

Y n (x) := 1 + Yj,nx- j , (2-46b) 



r w (a:) = ^ ^ r ^^ + T -i H T ) ' ( 2 " 46c ) 
k=i 

where the Tk are all diagonal and (in our case) traceless. In the present case we easily find T n (x) from 
the tracelessness and the exponential part of the asymptotic form of the fundamental system defined in 
eq. (2-31b) (cf. refs. || Q for similar formulas in the case of quartic even potentials). That is, 

T K = ^dmg(u K , -u K ) = ^a 3 , K>0, T_i = na 3 (2-47a) 

T n (x) = - (^V{x) - nln(*)) CT 3 • (2-47b) 

The Stokes matrices relate the basis of solutions with these asymptotic forms in contiguous sectors. The 
isomonodromic tau function is then defined by the deformation equation. (See the introduction to p2| 
for the general context and motivations behind the theory of isomonodromic tau functions.) 



dln(r 7M ,„) := res Tr (Y' 1 (x)Y^(x)dT (x)) = £ -L res Tr (Y- 1 (x)Y^(x) a^' 1 ) du K , (2-48) 
where 

d+i „ 



k=i 

We now state and prove the main theorem relating r/M,n( >£ 'i V) an( A Z n (>c, V). 

Theorem 2.4 The partition junction Z n (V) is related to the isomonodromic tau function by the formula 

TlM, n (^ V) = Zn(x, V) Ud+1 ^T , (2-50) 

up to a multiplicative factor that depends only on x. 
Proof. We will prove that 

o 

n 

dlnZ„ = d hi TiM,n ~ , ; ; tt du d +i 

(d + l)u d+ i 

We start by observing that (suppressing the dependence on x and n for brevity and denoting the derivative 
with respect to £ by a prime) 

V n = I-'*" 1 =CY'Y- 1 C- 1 +CYT'Y- 1 C~ 1 (2-51a) 
Tr(X>„ 2 ) = Tr ((y'y- 1 ) 2 ) +Tr Tt' 2 ) +2Tr(y- 1 y'T') . (2-51b) 

Consider the coefficient of x K in Tr(2?„ 2 ), which is: 

1 A d ^^ 1 d Thmj^ 

-} Uj +1 u k +i-j-2nu K+ 2-2} uj +k+2 Jt\ — ln(£„) — - res x K x Tr(T> n 2 ) = 

2 * — ' * — ' OU J x=oo 
3=0 J=\ J 
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=0 because Y'Y^ 1 =0(x- 2 ) 

'1 N 2 



= - res x- K - lf Tr ((Y'Y^ 1 ) 2 ) -2 res x~ K - 1 ±V'(x) - - - 2 res aT^Tr (y- x y' T') = 

a;— oo 2— oo \ 2 X / x—oo 

K d+1 

= 2 zJ J _ 2nu K+2 - 22, M ^ Tr o'as' 7- ^ -2 ) = 

j=o j=i ~ 

1 K d-Jf-1 g 

= 2 X! u .m u K+i-./ - 2nu K+2 - 2^ u J+K+2 J - — \n{TiM, n ) (2-52) 
J=0 j=i Uj 

We therefore have the relations: 

d-K-l g d-K-l q 

V" u J+K+2 J- — ln(Z„) = V" u J+K+2 J- — ln(r/ M ,™) , # = 0, ...d-2. (2-53) 
j=! j~i " UJ 



The vector fields appearing in eq. (2-53) are just representations of the negative Virasoro generators on 
the space of functions of the deformation parameters {uk}- 



d+l-K „ 

a 



J=l - JJ 



Eq. (2-53) can therefore be rewritten as 

V_ A - In (Z n (V)) = V_ A - In (r /M ,„(V)) , 2 < K < . . . d . (2-55) 

Since the vector fields Vo, ...,V_d are linearly independent 

V A ... A V_ d = {d+ 1)! {u d+1 ) d+1 d Ul A ... A 9 U , +1 , (2-56) 

and Z n (V), tim,u depend on the d+ 1 parameters iti, . . . , Ud+x, we may relate them up to a multiplicative 
constant if we determine the effect of the two remaining operators. Now Vo and V_i generate dilations 
and translations respectively 

Y V(x) = xd x V (x) , V_iK(a;) - d x V(x) . (2-57) 
Moreover it is clear from its definition (eq. (|2-§| )) that the partition function must satisfy 

V_iZ„ = , V Z n = ~n 2 Z n . (2-58) 



On the other hand the definition ( 2-4q ) of the isomonodromic tau function implies that it is invariant 
under translations and dilations of x so that 

V r 7M , n =V_ 1 r /w ,, = 0. (2-59) 

It follows that 

2 

71 

In (Z n ) = ln(T 7A f.n) - -7— ln(wrf+i) + const, (2-60) 
d+l 

which gives the proof of the theorem up to a constant multiplicative factor. Q.E.D. 

Remark 2.4 Thm. 2.4 was proved in the special case in which V[x) = g\x 2 + g2X 4 in | po| . 

Remark 2.5 For any fixed n, the set of >c for which orthogonal polynomials exist up to degree n is an open set 
because it is the complement of the vanishing locus of n homogeneous polynomials 

A m (x,y) / , Vm < n. (2-61) 

Thus we could also study the deformations with respect to parameters h. However, these are not isomonodromic 
deform ations because they change the data of the Riemann-Hilbert problem that is naturally associated to eq. 
(2-44). A similar approach in the case of even potentials can be found in 
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A Appendix 



In this appendix we give proofs of Theorems |2.1| and |2.2| . First, we consider the deformation equations, 
which determine the variation of the orthogonal polynomials when the potential V(x) is varied. 
Let us first expand the variation of ip n (x) with respect to uk in the basis {ip n }'- 

3 

i> n (*) = J2 V^nmTpM ■ (A-l) 



" m=0 



The semi-infinite matrix Uk is antisymmetric (as seen from differentiating eq. (2-11)): 

(U K )nm = -{U K )mn ■ (A-2) 

Taking the derivative of eq. ( |2-10| ) with respect to uk and using dp n /duK = 0(x n ^ 1 ) and eq. ( |A-2j )) we 
obtain 

tfe = ~((g*) + -(Q*)_) ■ ( A -3) 

In particular the semi-infinite matrix Uk is of finite band size: 

(U K ) nm = if \n - m\ > K , (A-4) 
and therefore eq. ( A-l ) contains only a finite sum. The diagonal part yields (from eq. ( 2- 10| ) ) : 

n 1 

i=0 

The semi-infinite matrix Uk can be used to form a differential system in the same way as for the d x 
equations, by expressing any ip m in terms of ipn and ipn-i by means of eq. ( 2-12 ). We then obtain an 
equation 

d f^-^\-U KAx) (^-^) , (A-6) 



du K \ 1pn(x) J V i>n(%) 

where Uk,ti(x) is a 2 x 2 matrix with polynomial coefficients of degree < K. In particular it is easy to 
compute U\^ n {x): 

U^) = -\(^: X ?\ V (A-7) 



2 I - 2 7n x - (3 n 



The general form of the matrices Uk,u is given in Thm. 2.2 which we now prove. 
Proof of Thm. |2.2| For any semi-infinite matrix A, let A\ denote the / th diagonal above (or below, if 
I < 0) the main diagonal. Recall that A + (resp. A_) denotes the strictly upper (resp. strictly lower) 
triangular part of A. 

Since the matrix Q has only three non-vanishing diagonals, we have: 

(Q K+1 ) + - (Q K Q)+ = {Q K Qi)+ + (Q K Qo)+ + (Q K Q-i)+ 

= ((Q K ) + + (Q K ) ) Qx + (Q K ) + Q + ((Q K ) + (Q K ) 1 ) Q-i 

= (Q K ) + Q + (Q K ) Qi-(Q A ') 1 Q-i • (A-8) 

This, combined with a similar calculation for (Q K+1 ) _ implies: 

Qt +1 - Q K+1 - (Q+ - Q K )Q + (Qo + Q^)Qi - (Qf + Q^Q-i , (A-9) 
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which, componentwise, reads (using eq. (A-3)): 

d d 

- 2{K + 1)- tpn{x) = -2K- X1p n (x) + (Q K ) rir , (jn+ltpn+1 ~ Jnlpn-l) 

ouk+i ' '" • 



Using eq. (2-12): 



du 



K 



(A-10) 



2(K+1)- 4> n (x) 

OUK+l 



■ ) 

-2{K + 1)- 7(>n-l(x) 

OUK+l 



= -2K—xMx) + (Q ) nn (0* - PnWn - 27nVn-0 



K 



+ {ln{Q K ) n>n _ 1 -ln + X (Q K ) n , n+1 )^n 



(A-lla) 



d 



-2K-£—ml> n -. 1 (x) + (Q Ar )„_ li „_ 1 (27n^« - (x - P n -l)lj) n -l) 



+ (7»-l (Q K ) n - hn - 2 ~ In {Q K ) n;n -i) ^n-X 
which can be further simplified using 

n,n+l n/n n,n—l 

Q n -l,n-l = ln-lQn-l,n-2 + Pn-lQn-l,n-l + lnQn,n— 1 • 

We thus get a recursion formula for Uk ,n (x) : 

-2{K + l)U K +i, n {x) + 2KxU K , n {x) = 



(xQ K -Q K +^) 1 1 

V ^ / n — l,n— 1 

(xQ* - Q^ 1 ) 



27n 



^rm ^n,n-l 



from which the result follows by induction on K and using eq. ( J A- 7| ) . Q.E.D. 

As a corollary to Thm. |2.2| we can complete the proof of Thm. [2.l[ 
Proof of Thm. 2.1. Since V(x) is a polynomial of degree d+ 1, eq. (2-20) implies that: 



(A-llb) 



(A-12a) 
(A-12b) 



(A-13) 



P = u k+i K U K , 



K=0 



and therefore 



K=l 



(A-14) 



(A-15) 



Using Thm. 2.2 and eq. (2-21), one immediately finds the expression (2-25b) for V n (x). Q.E.D. 



B Appendix 

We give here a simple proof of the fact that the partition function Z n (V) is a KP tau function in the 
sense of [3(J; that is, the determinant of a projection operator acting on elements of a suitably defined 
Grassmannian, evaluated along the orbits of an Abelian group. The interest in rederiving this well-known 
result is that the method used, which just involves a simple formal matrix calculation, naturally leads to 
a new interpretation of the partition function as the Fredholm determinant of a simple integral operator. 
Let 

H = H n +H+, (B-l) 
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where H n = spaxi{z 3 }o<j<n-i and H.+ is a suitably denned completion (whose details are irrelevant for 
this sort of formal computation) of the space span-f^ }.,■>„, with z£C viewed as a variable defined along 
our cycle xT. (Note that we could have included the negative powers {z°}j<o in the definition of H, but 
for the case at hand, these contribute nothing to the computation of the tau-function and it is simplest 
to suppress them from the start.) Denote by wq the Gl(TL) element defined by 



wo 
w 



H 

Sj>0 a 3 z:1 



H 

T,j>o a oPj( z ) 



(B-2) 



where {pj(z)} are the monic orthogonal polynomials associated to the moment functional C^r for an 
initial value of the parameters {uk}- Within the natural basis {z J }j£f$ the group element wo is represented 
by the upper triangular matrix 



(B-3) 



n 


+ 




' w n 


w n+ ' 


n 





W+ 


+ 



Wo := 

whose j-th column consists of the expansion coefficients in 

3 



P.i 



(z) : ^W kj Z k ■ 



(B-4) 



fe=0 



(Note that the blocks W n , W + are upper triangular with diagonal elements all equal to 1.) 

We consider the Grassmannian Gr + (H) of subspaces of Tt that are "commensurable" with Ti + C H 
in the sense of [E50l, and the commuting flows induced on it by the action of the Abelian group 



7(t) 



c Gl{H) 



(B-5) 



(where t = (ti, t2, ■ ■ .)), acting by multiplication: 

r+xH — > n 

(7(t),/) ^ 7(t)/ 



(B-6) 



Within the standard basis, the flows on Ti. are represented by semi-infinite lower triangular matrices 



7(t) 



" A 





B 


C ' 



(B-7) 



where the diagonal blocks are again lower triangular with l's along the diagonal. The initial element 
wo € Gr + {TL) is taken as 

wo :=spa3x{pj{z)}j> n . (B-8) 

In homogeneous coordinates, relative to the standard basis, this is represented by the (n + oo) x oo 
rectangular matrix 

Wo := Wo ° . (B-9) 

The induced flow on Gr+(H) takes wo to the element with homogeneous coordinates 

W(t) = 7(t)FT . (B-10) 

The corresponding KP tau function Tw (t) is given by the Pliicker coordinate corresponding to the 
subspace 

r Wo {t) = det {B W n+ + CW+) . (B-ll) 
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Since all diagonal blocks have unit determinants, this equals 



TWo(t) = det loo + BW n+ W+- L C- L = det 1 B + W^W+^C^B 



(B-12) 



where the second equality, expressing rvi/ (t) as a finite n x n determinant, follows from the Weinstein- 
Aronszajn identity. 

The partition function Z n in ( |2-§| ) is also expressible as an n x n determinant of the matrix of moments 

Wl 

9% := C„ T {x l+r ) , 



which in our notation may be expressed at the parameter values 

,(°) 



as 



where 



uj =uy + tj, J = 1, . . . , d+ 1 

t := (ii,i2, . . . ,id+i,0,0, . . .) , 
Tl = W^ 1 H W - 1 



is the initial value of 971 (at t = 0) and 

H := diag(/i , hi, . . .) 



Ho,n 
^o,+ 



is the initial value of the matrix of normalization factors. 
Let 

hi '■ 7~L n ^ % , 7Tn : H — > H n , 

be the natural injection and projection represented by the matrices: 



7r n := [1„ 0] 



Then, by eq. (2- 



±Z n = det(ir n VJli n ) = det ^„W 4 H W Q - S" 
= det(i? 0i „)det(l„ + Fy„ + ^+- 1 C- 1 B) , 
where the second line follows from the fact that 

det (A) = det(C) = det(W n ) = det(W+) = 1 . 

Therefore 



rw (t) = 



Z n {t) 



Q.E.D. 



(B-13) 

(B-14) 
(B-15) 

(B-16a) 
(B-16b) 

(B-17) 

(B-18) 
(B-19) 

(B-20a) 
(B-20b) 

(B-21) 
(B-22) 



Zn(0) 

We may now use this calculation to deduce a representation of Z n (t) as a Fredholm determinant. First 
note that from the definitions of the matrices A, B, C, W n +, W+, we may view these as representations 
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of the operators: 



B : Hn 

C : H+ 
( ' : Ej h 
W+ : H+ 
>'»'• : Ej>„ 

: n+ 

W n+ ■ T,j>n a 3 z3 



E"=o^(7(t)^)>„ 
Ej>„ a j7(t)z J 

■ n+ 

Ej>n a ](Pj( Z ))> n 

* Ej> n aj(pj(z)) <n 



(B-23a) 
(B-23b) 
(B-23c) 
(B-23d) 



where (f{z))> n and (f(z))<n denote the parts of the Taylor series expansion of / G TL lying in H.+ and 
H n respectively. It follows that the operator : 



K : H+ — > H+ 



represented above by the matrix appearing in eq. (B-12): 
is just the rank-n Fredholm integral operator 



i " <0> 

Kf(z)= I dwe-Z d *=i-B- wK K(z,w)f(w) , f G H+ 

JitV 

n-1 

K(z,w) = E - (e^ t ^ K / K p j (z)) >n P j (w) 
7w„(t) = det (l-^) • 



with kernel 



and hence, by eq. ( B-12 ) 



d+1 . K I js 

L K=1 t K w /K 



(B-24) 
(B-25) 

(B-26) 
(B-27) 
(B-28) 
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